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Abstract 

In this paper, we study the physics of mesoscopic systems with noninteract- 
ing, but fixed number of electrons. From a technical point of view, this means 
a discussion of the differences between the canonical and the grand canonical 
ensemble (fixed versus fluctuating number of particles). Such a discussion is 
not trivial since the grand canonical ensemble is the most convenient basis 
for the statistics of identical particles and one has to spend labour in order 
to retrieve the canonical ensemble. Specifically, we are considering ensembles 
of mesoscopic systems with disorder, either by atomic defects or by fluctua- 
tions in their geometric definitions and we discuss various forms of disorder 
averages. 
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I. INTRODUCTION 



In the following, we will investigate some properties of noninteracting electrons in meso- 
scopic systems. We take into account that irregularities in the preparation introduce disorder 
on an atomic scale and also, in the geometrie definition of the samples. We will explore the 
consequences of strict conservation of particle number and we will find that some details 
in the disorder averaging procedure may become important. Such a peculiar behavior has 
been emphasized by various authors; we wish to mention here Shklovskii |T| as well as Imry 
0. From a technical point of view, our paper is meant to discuss the differences between 
the canonical and the grand canonical statistical ensemble. 

We will be concerned first with (A) thermodynamic properties on the average ||, sec- 
ondly with (B) kinetic behavior and thirdly with (C) stochastic properties of thermody- 
namic quantities [3]. Clearly, in the first and third case, it is suffice to know the positions 
of the single particle levels €\. We will argue that the position of the levels is also most 
decisive in the second case. 

(A) According to standard arguments, thermodynamic properties may be derived from 
the grand canonical potential f2(T, /i) which is given by 



where s = (2S + 1) is the spin degeneracy. (In previous publications, we used to start with 
the density of states per spin and unit volume Af(E) = V(E)/V, following a convention of 
previous decades, where one has been interested in bulk properties.) Then, the expression 
for the grand canonical potential can be written as 




(1) 



A 



We introduce the single particle density of states per spin 
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(B) Following Ref. U, we will also study the photoabsorption of mesoscopic particles. 
The cross-section for this process can be written in the form 

47r 2 e 2 
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where ri\ is the population of the single particle state | A >; and eR\y = e < A | R \ A' > 
the matrix element of the electronic dipole moment. It seems (see Refs. [4] - [6]), that there 
is not much structure in | R\y | 2 so that one may replace it by an average value 

47r 2 e 2 l9 

(which also includes an average with respect to the disorder), without distorting the general 
structure of the result to be obtained. Assuming that the population is equal to the Fermi 
function, we obtain 

a = a Q h 2 uo s J dEdE'V(E)V(E') f(E - /i) - f(E' - //) 5(ftw - E' + E) (6) 

Due to the discreteness of the levels in mesoscopic samples, the tu-dependence of the cross 
section <j(uj) consists of 5-spikes for each individual sample and it is only the average with 
respect to the ensemble of samples, that is, the disorder averaged cross section < criu) >, 
which may be expected to be a continuous function of the frequency. 

As a rule, the density of states is a fluctuating quantity which consists of a leading part 
V°(E) and a remainder 

V{E) =V°(E) +V\E) (7a) 

where 

V°{E) =< V{E) > (7b) 

is equal to the average with respect to the disorder. Note that in the literature, the following 
abbreviation occurs quite frequently 

TZ 1 {E) = V°{E) (7c) 



(one level correlation functions); and if systematics is required, we will resort to that termi- 
nology. 

We mention in passing that in some problems where no disorder appears explicitly, 
the average can be taken with respect to an energy window or it may represent a kind of 
integration with respect to quantum numbers. 

As a rule, T>° is insensitive to genuine quantum phenomena which may result e.g. from 
an applied magnetic field or from a magnetic flux threading a mesoscopic ring. 

If this is the case, the disorder averaged grand canonical potential 

n° =:< n>=s J dEV°(E)(-T)£n[l + e ~ {E ~^ /T ^ (8) 

has lost any sensitive dependence on external parameters. Indeed, this turns out to be true 
as far as the response (persistent currents) to magnetic flux of mesoscopic rings is concerned. 
On the other hand, it is found that there is a measurable effect in a theory where the number 
of electrons is kept fixed [2,7]. 

(C) The fact that the sensitive part of Q is, on the average, zero does not mean that Q 
is zero for each realization. Rather, it means that its contributions fluctuate in magnitude 
and in sign; loosely speaking, one may say that plus and minus signs occur in about equal 
numbers. The simplest quantity which can be used as a measure of the stochastic (sample 
to sample) fluctuations, is the mean square value. In previous publications (see references 
in [3]), it has been considered to be sufficient for non- interacting electrons, to calculate 
the mean square < Q 2 >, which means stochastic fluctuations at fluctuating numbers of 
electrons. We will show how stochastic fluctuations can be calculated for fixed number of 
electrons. A preliminary estimate of the resulting expression indicates that significantly 
larger values are then obtained as compared with the standard procedure. 

For an initiation to a discussion of the photoabsorption problem, we calculate the disorder 
averaged cross section in a crude approximation < a > uc where in Eq. (|6|) the averaged 
product of the density of states is replaced by the product of the averages (uc: uncorrected; 
that is, disregarding correlations). If Ttw « [i, we may even put V°(E) where 



P V) == V F 



(9) 



is the average density of states at a (typical) Fermi level Then, we obtain 



< a > uc = h 3 u 2 sV 2 F a 



(10) 



The paper is organized as follows. Thermodynamic properties (A) of mesoscopic systems 
with fixed number of electrons are investigated in Section II, III, and IV on the basis of (i) 
Legendre Transformation; (ii) Coulomb blockade; and (iii) pinning of the Fermi level to 
a single particle state. Section V is devoted to the problem of photoabsorption, whereas, 
stochastic fluctuations are studied in Section VI. The paper concludes with a discussion 
in Section VII. Appendix A presents a calculation of correlation functions in the cooperon- 
diffuson approximation for persistent currents and Appendix B contains a collection of results 
of the random matrix theory. 

II. CANONICAL ENSEMBLE BY LEGENDRE TRANSFORMATION 

The mathematical conveniences provided by the grand canonical ensemble (fluctuating 
number of particles) are evident and there is no need to consult text books. Certainly, there 
are differences in the results obtained for the above mentioned ensemble and the canonical 
ensemble (fixed number of particles) but they are said to become relatively small in the 
thermodynamic limit. 

However, there are phenomena in mesoscopic systems which depend only on a relatively 
small number of electrons. This suggests that one should work entirely within the canonical 
ensemble; in fact, this retriction should be taken very serious since in all processes the 
particle number is conserved. 

At low temperatures (at T = 0, to be precise) where thermal fluctuations are negligeable, 
the grand canonical potential f2(T = 0,/i) and the free energy F(T = 0,N) are directly 
related by a Legendre transformation. However, this is true only before the averages with 
respect to the disorder have been taken; therefore, the transformation has to be done for 
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each realization of the disorder separately. This program can be carried through comparably 
easy, on the basis of the decomposition (7), provided that the fluctuating part T> l (E) in the 
density of states is small. We will find that the first nontrivial correction is of second order 
in D 1 . In what follows, we will relax the condition T = and we will perform the Legendre 
transformation from £l(T, fi) to F(T, N) on the basis of the above mentioned expansion [7]. 
In this context, it is convenient to formulate the Legendre transformation as follows 



F(T, N) = max^|fi(T, //) + /jln} 



(11) 



Clearly, a necessary condition for a maximum is 

+ jV = ^ /i = /i(T, N) 



(12) 



Corresponding to the decomposition (7) of the density of states we have Q = Q° + Q 1 , where 
Q° =< Q > is given in Eq. flD and where 



Q 1 = s rfEP 1 (E)(-T)£n 



1 + e 



(13) 



Let us put 



(14) 



where /i° is the zero order approximation such that 



N 



(15a) 



/'o 



and let us assume that we have found the leading order Legendre transform 



F°(T, N) = fi°(/i°) + fi°N 



(15b) 



According to what we have said in the Introduction, F° does not depend on the quantum 
effects we are interested in; therefore, its detailed form is not important here. In the next 
order, we have 



on 

djj, 



/i° + /i 1 



on 1 

djji, 



N 



(16a) 



Making use of the relation (cf. Eq. ( 

d 2 n° 



dN 



= sV F (16b) 

o 



9/i 
we obtain 

sVpfi 1 + s J dEV 1 (E)f(E - n°) = 
f=±- f dEV 1 (E)f(E — fj?) (17) 

We conclude that through second order 

F = F° + Q 1 {fi°) + AF (18) 

where 

^^WW + fliT"} (19) 

and where the derivatives are taken at fi = fi°. Consequently, 

1 tfW, lN2 



2 9// 
1 1 



2 



s / dEV 1 (E)f(E - fi°) (20) 



2s£> F 

For convenience, we introduce here 

A = 1/sVp (21) 

which is the mean level separation and 

SN = s J dEV\E)f{E - n°) (22) 

which is the disorder induced fluctuation in the particle number calculated for the grand 
canonical ensemble. Thus, we may write 

AF = ^A(5N) 2 (23) 



7 



We recognize that the calculation of the disorder average < AF > requires the knowledge 
of the two-level correlation function 

TZ 2 (E h E 2 ) =< V(E 1 )V(E 2 ) > (24a) 

which is related to the cumulant as follows 

y 2 (EiE 2 ) = -n 2 (E 1 ,E 2 )+n 1 (E 1 )TZ 1 (E 2 ) 

= - < V 1 {E 1 )V\E 2 ) > (24b) 

For orientation, we mention some properties of this correlator. We expect that for large 
separation of the energies the correlations vanish. Hence, 

J dE y 2 {E,E') — (25) 

In general, we expect only a weak dependence on the absolute values of the energies; hence, 
for energies near the Fermi level, the correlator may be assumed to depend only on the 
energy difference 

y 2 (E,E')=Y 2 (E-E f ) 

= Y 2 {\E-E'\) (26) 

where the second line follows from the symmetry E <-> E'. 
Making use of Eq. (25) and then of Eq. (^), we arrive at 

< (SN) 2 > = s 2 J dEdE'f(E - fi°)f(-E' + fi°)Y 2 (E - E') 

1 / l + e u ' 2T 



= S ' 2T J dU -^{^ l + e -u/2T ) Y ^) (27) 

In the limit T — > 0, the expression simplifies to 

r° 

< (SN) 2 >=-s 2 duuY 2 (u) 

= s 2 duuY 2 (u) (28) 



o 



In retrospect, one comment is in order. We should recognize that the two energy integrations 
in the expression for < (5N) 2 > collect contributions from a large range of order Fermi 
energy. Therefore, it is not obvious that the translational invariant form (26), for the 
correlator, is admissible. For justification, note that firstly, the application of the sum rule 
(25) allows to introduce substantial restrictions in one energy integration. Still, the one 
energy integration in Eq. (28) may diverge. There, we should keep in mind that, we are 
interested here in quantum effect which occur at low energies. Therefore, it is possible to 
subtract irrelevant parts and to arrive at a convergent expression one is interested in. 

For a demonstration, see the discussion on persistent currents, where a metallic ring is 
threaded by a magnetic flux <f> (see appendix A). There, it is found that the phase sensitive 
part of < AF > is periodic in <fi with period where O — 2Trhc/e is the flux quantum. 
It is also found that the Fourier components of < AF >, that is, the phase sensitive parts 
of < AF >, are of the order of the mean level spacing A. 



Here, we discuss the consequences of the fact that the electrons carry the charge (— e) 
and that the capacitance C of a mesoscopic sample is small. (In case of a spherical geometry, 
C is proportional to the diameter.) 

Suppose that there is an excess charge (—e)(N — iVj) where N is the number of the 
electrons and iVj the effective number of the positive ions. Then, the electric potential is 
non-zero and equal to 



III. CANONICAL ENSEMBLE BY COULOMB BLOCKADE 



<j> = C(-e)(N-Ni) 



(29) 



Accordingly, the single particle levels undergo a change 



e A -> e A - e0 



(30a) 



which can alternatively be expressed as a change in the chemical potential 
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ji — > n + e0 (30b) 

For a discussion of the thermodynamics of the system, it is most convenient to use 
instead of N as an independent variable. In this case, the charging energy C0 2 /2 appears 
with a minus sign in the expression for the thermodynamic potential; thus 

fi(T, /U ,0) = fi(T, / x + e0)-ic0 2 (31) 

where the "free" thermodynamic potential Q(T,fi) is given by Eq. (3). 
In thermal equilibrium, f2(T, /i,0) is maximal with respect to 0; hence 

ft c (T, At) = min fi(T, /i, 0) (32) 

and the ensuing argumentation is similar to the one below Eq. (11). Clearly, a necessary 
condition for a minimum is dVL(T, /i, </>)/d<f> = 0; that is 

- e (N(T, h, 0) - ivV) - C0 = (33) 

where the contribution of the positive ions has been subtracted "by hand" . 

Following the decomposition (7) of the density of states, we have f2° + il 1 as well as 
N = N° + N 1 . Let us now assume that 

iy°(T,/i,0 = O) = iV, (34) 

Note that this assumption corresponds to what has been expressed by Eq. (15a) where the 
"typical" Fermi level /x° is defined; for convenience, we have omitted the superscript zero 
and we have put /i — /i° in (34). 

As previously, we consider the fluctuating contribution N 1 (which is proportional to V 1 ) 
as well as to be relatively small. As a consequence, Eq. (31) can be expanded as follows 

n(T, /i, 0) = n° - N 1 ■ (e0) - - ^(e0) 2 - -C0 2 (35) 

where the arguments of the thermodynamic quantities on the right side are (T, //). Calcu- 
lating the minimum of (35), we find 
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tt c = tt° + AQ C (36a) 

where [8] 

^•^^W^'M 2 < 36b > 

In the context of the Thomas- Fermi theory, one may consider the replacement C/e 2 — > 
C/e 2 + sV F as a manifestation of screening. 

In the limit where the charging energy Eq is much larger than the mean level distance, 
that is 

we arrive at 

Afi c (T, / i) = iA(iV 1 (T, / i)) 2 (38) 

which agrees exactly with AF of Eq. (23) if the identity N 1 = 5N - see Eq. (22) - is taken 
into account. 

The physics which has lead us to Eq. (38) has been called in Ref. [7] global charge 
neutrality. We wish to add that it is quite legitimate - if not obvious - to consider the 
mesoscopic sample connected to an electron reservoir. Of course, differences in the work 
function between reservoir and sample - which occur in a real situation - may require a 
proper redefinition of the typical Fermi level n°. 

IV. PINNING OF THE FERMI ENERGY TO A SINGLE PARTICLE LEVEL 

The theory which we will put down below is based on a statistical concept of the filling 
factors of a Fermi system in the presence of disorder. At T = 0, that is, in the ground 
state, all single particle levels are fully occupied which satisfy t\ < ep] the remaining levels 
(ca > €f) are unoccupied. Of course, the Fermi energy ep has to be chosen such that the 
particle number N which follows from this choice, satisfies the specifications. 
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At that point one may ask: what are the specifications? In a mesoscopic sample of size, 
say 10 3 atomic distances, is there the number of electrons exactly 10 9 , or perhaps 10 9 + 10 2 
(10°, 10 1 , 10 3 , ...)? Physical considerations suggest that the result should not depend on 
such details. 

Therefore, we propose to average with respect to the particle number in a reasonable 
window. A convenient way to realize such an ensemble will be to select Fermi-energies at 
random such that 



where e\ may be any level in a reasonable range (which will be specified below); and where 
5 is a positive infinitesimal [9]. The above choice implies that the levels e A < e A only 
are occupied; and that the remaining ones are empty. We will call this construction the 
Fermi-level pinning ensemble (FLPE). 

At that point, one may ask: what is the difference then of this ensemble and the grand 
canonical ensemble? The answer is implicitly given by arguments found in Ref. [1]. Ac- 
cordingly, the grand canonical ensemble (GCE at T = 0) may be characterized by a random 
selection of the chemical potential. This implies that the transition from the last occupied 
level to the first unoccupied one will preferably take place when the energetic distance is 
large. In contrast to it, the ensemble (FLPE) introduced above guarantees that such a 
transition occurs between any pair of levels with equal probability. 

Again, we start from the expression (1) for the grand canonical ensemble; however, we 
wish to add that only the limit T — » (ground state) is, strictly speaking, consistent with 
the argumentation above. Inserting relation (39), we obtain the following relation for the 
"pinned" grand canonical potential 



A 

As a rule, the average properties of a Fermi-system do not depend on the absolute values 
of the single particle energies. Therefore, we remove the arbitrariness in the selection of the 



€f = e\ + 5 



(39) 




(40) 
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pinning level by sampling with a weight function P(e\) which is centered at the typical 
Fermi energy, that is, = /jP (see Eq. (15)) and which has a support much smaller than /i 
but much larger than the mean level spacing A. Thus, we obtain for the average potential 

< V P >= < ^ p( ^ > < J2 P (ex)(-T)en[l + e-^- s » T ] > (41) 

AA 

Note that we have taken the average separately for numerator and denominator. One may 
justify this procedure in view of the assumed properties of the weight function P(e^). Nev- 
ertheless, this procedure may give rise to delicate questions since the properties we are 
interested in may depend on one electron only. 

Inserting the definition (2) of the density of states, and observing the relations (7) and 
(24), we obtain 



< Qp > 



J dEdE P{E) TZ 2 (E, E)(-T)ln 1 + e -(^~ 5 )/ T ] 



/ dEKi{E)P{E) 

= < Attp > (42) 



where Q° is defined by Eq. (8) and where 

< AQ P >= - d z n ° { E )P{ E) J dEdEP(E)y 2 (E,E)(-T)£n[l + e-( E - E - 5 » T ] (43a) 

By and large, the E integration above collects contributions from a large range of order Fermi 
energy. Therefore, we may not insert at once the translational invariant form (26) for the 
correlator. On the other hand, the E integration covers only the small range | E — /x° |~ A 
and there, we may assume an appropriate form of translational invariance. Thus, we may 
write 

< AQ P > = VfJ ~ S Ep{E) J dEdEP(E)y 2 (E,V° - S)(-T)£n[l + e'^^] 

= W F J dE ^ E ^° + S)(-T)en[l + e~ E -^' T ] (43b) 

Note that for T -> 0, we have (-T)£n[l + e -^-f°)/T] _^ e(/2° _ E ^o _ E y Therefore; the 
small quantity 8 is irrelevant. Furthermore, since we are interested only in quantum effects, 
which occur at low energies ("phase sensitive contributions"), we may now make use of the 
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translation invariant form (26). Introducing the integration variable u = E — /x° we arrive 
at 



< AQ P >= -As 2 / duY 2 (u)(-T)£n 



1 + e 



-u/T 



(44) 



Surprisingly, we recognize that 



< Afi P (T = 0) >= -As 2 j duuY 2 {u) (45) 

is just twice as large as < AF(T = 0) > as given by Eqs. (|3|) and (^8|). This apparent 
discrepancy can be understood as follows. We recall that essentially one is interested in the 
dependence of Q = Q(y) on an external parameter y. Such a dependence is given when the 
single particle energies e\ = €\(y) also depend on this parameter. Thus 



dy ^— ' dy 

A 



(46) 



Clearly, the y-dependence will also show up in the correlation functions. Therefore, we 
define 



where 



ME, y;E,y) = -< V\E, y)V\E, y) > 



V(E,y) = -J2S(E-e x (y)) 



(47) 



(48) 



and V 1 = T> — D°. (We repeat that T>° =< T> > is not expected to depend on y in an 
essential way.) For energies not far from the Fermi energy, the correlator will be stationary 
with respect to energy 



ME,y;E,y)=Y 2 (E-E ] y,y) 
= Y 2 (E-E;y,y) 



(49) 



and thus, we may generalize Eq. ([44]) as follows 



< An P (y,y) >= -As 2 / duY 2 (u\ y, y)(—T)£n 



1 + e 



-u/T 



(50) 
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According to what has been said above, the quantity we are interested in is 



< X(y) > 



d 



< n P (y,y) > 



(51) 



dy 



y = y 



Now, if Qp were symmetric in y,y f , we would be able to write 



< X(y) > 



ld_ 
2dy~ 



< tt P (y,y) > 



(52) 



and to remove thus the discrepancy, at least as far as physically measurable quantities were 
concerned. 

Comparing the structure of < AQ P > and of < 5N(y)5N(y) >, one recognizes that such 
a symmetry exists at T = 0. For T^O, one finds such a symmetry for Y 2 (E — E; y, y) in the 
case of a metal ring threaded by a flux (persistent currents; see appendix A). Of course, one 
should always keep in mind that both procedures (Legendre transformation and pinning of 
the Fermi-level) have a sound basis only for T = 0. 

In conclusion, we have to point out one peculiarity of the Fermi-level pinning ensemble. 
We should not expect any change in the final result if the selection ef — ^\{y) + S of Eq. 
(39) is replaced by e F = e^(y) + 5. This expectation is supported by the fact that the energy 
levels as a function of y do not cross in disordered ("chaotic") systems. Therefore, the set 
of levels defined by P(e\(y)) ~ should be the same as P{ex{y)) provided that the weight 
function P is sufficiently broad on a range | e\{y) — £\(y) |~ A. However, one finds in the 
case of persistent currents such chances (see appendix A). Presently, we cannot understand 
this dependence on such a detail in the preparation of the Fermi-level pinning ensemble. 

For sake of completeness, we wish to mention in this context that Kamenev and Gefen 
[10] have discussed a gedanken experiment, where the system is prepared in a grand canonical 
state at y whence it is transferred adiabatically to the parameter value y (see the discussion 
in appendix A). 
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V. PHOTOABSORPTION WITH FERMI LEVEL PINNED 



In what follows, we will discuss the consequences of the ansatz (39) when inserted into 
relation ([]) for the photoabsorption. In a first step, we replace /i = ep = t\ + 5 in the 
argument of the Fermi functions. Next, we take an average with respect to the pinning 
levels by sampling with the weight function -P(e^). Thus, we obtain the following expression 
for the photoabsorption 

"= JimmW) J <tE«e d EV { W(E>)V { E)P { E) 
x f(E-E-5)- f(E' - E - S)\s(huj - E' + E) (53) 

Considering the disorder average < a >, we argue as previously, that it is possible to 
average numerator and denominator separately. In this process of averaging, the three-level 
correlator appears 

TZ 3 (E U E 2 , E 3 ) =< V(E 1 )V(E 2 )V(E 3 ) > (54) 

appears. Since we are interested only in energies close to the typical Fermi energy fi°, the 
correlator TZ 3 may be considered to be invariant with respect to a translation along the 
energy axis. Thus, we obtain 

<a > = ct -^s / dEdE"}Z 3 (E, E', ^-8) 
Df J 

x f{E-(j?) -f{E'-fi°)]5{hu-E' + E) (55) 

For sake of simplicity, we take the limit T — > where the Fermi function becomes a step 
function. In this case 

<a> =a s du1l 3 (u + /J,°,u + hw + - 5) (56) 

J-hw 

If we were to disregard correlations, 1Z 3 — > T> F , we would recover Eq. (10). 

It will now be necessary to study some general properties of this correlation function. 
Firstly, we define the three level cumulant as follows 
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y 3 (Ei, E 2 , E 3 ) = H 3 (Ei, E 2 , E 3 ) 

+n 1 (E 1 )y 2 (E 2 , e 3 ) 
+n 1 (E 2 )y 2 (E 1 , e 3 ) 
+n 1 (E 3 )y 2 (E 1 , e 2 ) 

-^1(^1)^1(^)^1(^3) (57) 

Again, we expect that for large separation of the energies the correlations vanish; hence 

J dE 3 y 3 {E u E 2 ,E 3 ) = (58) 

It is possible to find definite forms for the correlation function within the random matrix 
theory [11,12]. In our notation, we follow most closely the article of Bohigas [12]. 

However, there is one feature, the autocorrelation namely, which does not seem to have 



received sufficient attention in the past. For instance, consider Eq. ( |24a| ) and insert definition 
(§). One obtains 

K 2 (E, E')= K 2 (E, E') + 8(E - E')H X (E) 

TZ 2 (E, E')= -^J2 <6 ( E ~ E ^ E - E ^ > ( 59 ) 
Correspondly, Eq. (24b) is of the form 

y 2 (E, e')= -n 2 (E, e') + n^n^E') 

= y 2 (E : E') - 8(E - E')lZi(E) (60) 

such that 

j dE'y 2 {E,E') = K 1 {E) (61) 

Autocorrelations play an important role in photoabsorption. In thermodynamics, however, 
they are irrelevant as one may convince oneself by inspection of Eq. (28) and Eq. (44). 

For energies close to the center of the energy band, that is, for energies close to the Fermi 
level 
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we have 



£ =E-/2° 

(62) 

£' =E'-fi° " 



H^E-fi ) =i?i(5~0) =V F (63) 

Furthermore 

y 2 (S, J E?')=y 2 (5-5') 

F 2 (£ - £')= Y 2 (£ - £') - S(£ - 5')i2i(0) (64) 

where the following rules of notation have been introduced: Correlators with script letters: 
arbitrary energies Ef, correlators with roman letters: small energies £j. In this limit, the 
correlators are invariant with respect to translations. Correlators with (without) tilde: with 
(without) autocorrelations. 

For reasons of simplicity, the following realtions pertain to the Gaussian unitary ensemble 
(GUE) of the random matrix theory. There, one has 

Y 2 {£ - £') = i?2(o) s 2 (x - x ') (65a) 

where 

sin Tix 
s{x)= 

TIX 
£ 

Xi= (65b) 
We convince ourselves by explicit calculation that 

J d£' Y 2 {£-£') = i?i(0) (66) 

as it should be. 

A detailed inspection shows that the three level cumulant of Eq. (57) 

%{E^E 2 , E 3 ) = Y 3 (£ u £ 2 ,£ 3 ) (67a) 
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comprises the following (^-function contributions 



^3(^1, £<=, £3) — ^3(^1, £2, £3) 

-5(S 1 -S 2 )Y 2 (S 2 -S 3 ) 
-5(£ 2 - £ 3 )Y 2 (£ 3 - £ ± ) 
-5(£ 3 -£ 1 )Y 2 (£ 1 -£ 2 ) 
-26(£! - £2)6(82 - £ 3 )Ri(0) 

Using the above decomposition, we conclude that Eq. (58) assumes the form 



j d£ 3 Y 3 (£ 1 ,£ 2 ,£ 3 ) = 2Y 2 (£ 1 - £ 2 ) 



For the GUE specifically, we have 



Y 3 (£i, £ 2 , £ 3 ) = R 3 l (0)2s(x l - x 2 )s(x 2 - x 3 )s(x 3 - xi) 



For a control, we calculate 



dx 3 s(x 2 - x 3 )s(x 3 - xi) = s(x 1 - x 2 ) 

and thus, we find Eq. (68) confirmed. 

In accordance with the notation introduced above, we put 

TZ 3 (Ei, E 2 , E 3 ) = R 3 (£i, £ 2 , £ 3 ) 

and find the explicit form 

R 3 (£ 1 ,£2,£ 3 )=Y 3 (£ 1 ,£ 2 ,£ 3 ) 

-i?!(0){y 2 (£: 2 - £ 3 ) + Y 2 (£, - £ 3 ) + Y 2 (£, - £ 2 )] 
+5(£ 1 -£ 2 )[r 2 1 (0)-Y 2 (£ 2 -£ 3 ) 
+5(£ 2 -£ 3 )[r 2 1 (0)-Y 2 (£ 3 -£ 1 ) 
+5(£ 3 -£ 1 )[r 2 1 (0)-Y 2 (£ 1 -£ 2 ) 
+R 1 (0)5(£ 1 - £ 2 )5(£ 2 - £ 3 ) + R\{<S) 
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We insert this form in the relation (56) for photoabsorption. Since < a > is an even 
function of u, we take uj > 0. Firstly, we recognize that the contribution with a product 
of two ^-function in Eq. (71b) drops out. Next, we also recognize that the infinitesimal S 
is important for an exact definition of the contributions, which the remaining (^-functions 
collect from the end point of the integration interval. 

Using the specific forms of Eq. (65) and Eq. (69) and normalizing the cross-section 
< a > to the uncorrelated one of Eq. (10), we may put the result in the following form 

<a> 2 , , f x/2 
-six 



< a > uc x 



-s 2 (x) 



where 



f x/z x x 

J s{y+-)s{y--) 

+ 1 _ l s ^ x ) _ I f d y s 2 (y) + 1 (72a) 

XX X J _ x 



< 72b » 



For small and large frequencies, the above relation reduces to 



< a > 



(73) 



(2rr 2 /3)ftLu/sA fko « A 
< ° >uc [l + sA/fru fku»A 
By numerical integration, we have calculated the cross-section for intermediate frequencies. 
The result is shown in Fig. 1 

The same type of analysis can also be done for the Gaussian Orthogonal Ensemble (GOE) 
of the random matrix theory. This is outlined in appendix B. In this case 
<a> f (tt 2 /3) - (tt 4 /30)(^/sA) 2 ftw « A 



<a>™ { i + s A/hw hw»A 
For intermediate frequencies, the cross-section is found by numerical integration; for a graph- 
ical representation see Fig. I. 

The limiting forms (73) and (74) are in agreement with Shklovskii's [1] conclusions. 
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VI. STOCHASTIC FLUCTUATIONS 



Up till now, it seems to be an accepted procedure in mesoscopic physics to calculate 
stochastic fluctuations within the grand canonical ensemble. Thus, one calculates the con- 
nected correlator of the grand canonical potentials (see also Eq. (47)) 



< n(T,n,y)£l(T,n,y') > c =< Sltf > - < SI >< St > 



= —s 



J dEdEy 2 (E,y;E',y')T 2 (£n 



1 + e 



-(E-n)/T 



(e n [l + e-( E '-M T ]) (75) 



According to Eq. (46), the correlator < X(y)X(y') > c for the physical quantity X(y) can 
be obtained by operating with d 2 /dydy' on Eq. (75). 

As an illustration, such a calculation is presented in appendix A for the case of persistent 
currents. One obtains for the phase sensitive part of < QQ' > c by order of magnitude ~ E c , 
where E c = D/L 2 is the Thouless energy in the diffusive limit. Note that E c » A in 
experiments of the usual type. 

However, with the results and with the insight obtained in the discussions of the preceding 
sections, we may quite well ask what the stochastic fluctuations are when calculated for the 
Fermi level pinning ensemble. Considering Eq. (40), we write down in a first step 



n P (y,y)n P (y',y) = $> 2 (fn[l + e-^ +& ^ T ) [in 

AA' 



1 + e 



) (76) 



Next, we perform an average of in a reasonable energy range on the basis of the sampling 
function P(e\)- As previously, we claim that it is possible to take the disorder average for 
numerator and denominator separately. Thus, we obtain 



< n P (y,y)n P (y',y) > 

o2 



V F J dEP(E) 



dEdE'dEP(E) < V(E,y)V(E',y')V(E,y) > 



xT 2 Un 



1 + e 



-{E-E)/T 



in 



1+e -(E'-E)/T^ 



(77) 



As far as the energy integration is concerned, we wish to recall the argumentation in 
connection with Eq. (43). Accordingly, we write (typical Fermi-energy /i° — > /i) 
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< n P (y,y)n P (y',y) > = ^J dEdE ' < V{E,y)V{E',y')V{ii + S,y) > 

xT 2 (in [l + e -( B -")/ T ] ) (in [l + e^'"^] ) 



(78) 



Since < Q(fj,,y)Q(fj,,y') > can be obtained rather easily, it is advantageous to calculate 
the difference 

< n P (y,y)n P (y',y) > - < Q(fi,y)Q(fi,y') > 

= s 2 J dEdE'[^-TZ 3 (E, y- E', y'; fx + 5,y)- K 2 (E, y; E', y')] 
xT 2 (in (in [l + e~ (E -^ T ] ) (in [l + e^ E '-^ T ) (79) 

For a definition of TZ 2 and Tl 3 , see Eqs. (24) and (54), respectively. As far as the difference 
of the correlators in the curly brackets is concerned, we make use of Eq. (57) and obtain 



—ME, y, E', y'; n + 8, y)-K 2 (E, y; E', y') 
Up 

= -Ly 3 (E,y;E',y';ii + 5,y) 
Up 

-^-Ri{E)y 2 (E', y'- n + 6,y)- ^-(E')Y 2 (E, y^ + S, y) (80) 

Up Up 

As emphasized repeatedly, we are considering systems where R\(E) =< V(E,y) > does not 
depend on the external parameter y. 

We insert the last term of Eq. (80) in Eq. (79) and find 

„2 



-— J dEdE"Jl 1 (E')y 2 (E,y; fi - 5,y)T 2 (en 
= n (fi)<AQ P (y,y)> 



1 + e 



-(E-ix)/T 



In 



1 + e -(^M)/Tjj 



(81) 



The last line follows from Eqs. (8), (43) and (50). Similarly, the second last term of Eq. (8) 
contributes with < Att P (y', y) >. 

At this, we recall that ultimately, we are interested in the correlator of the physical 
quantity X(y). This means that we should calculate 



rj2 

< X{y)X{y\> P) y= — < n P ( y ,y)Q P (y',y) > 



(82) 
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(The subscript P indicates how the average is taken and y is the parameter value at which 
the Fermi level pinning ensemble has been prepared.) It is now important to realize that 
terms of the type (81) do not contribute to the correlator (82). 

For sake of transparent formulae, we will omit terms of the type (81); in the following 
equalities are meant to be relations of equivalence in the sense of Eq. (82). Therefore, we 
may write 

< n P (y,y)n P (y',y) > - < Q(fi,y)n(fi,y') > 



s*_ 



J dEdE'% (e, y; E', y'; (1 + 8,$) T 2 (in [l + e"^"^] ) (in [l + e^'"^] ) (83) 



Note that for T — > 0, the last factors in the integrand are 88'9(—8)9(—8') where 8,8' 
are energies measured from the Fermi level as defined in Eq. (62). Clearly, large values of 
e, e' are important in the integral (83), and a simple RMT theory fails if not a prescription 
how to extract phase sensitive contribution is supplied. As an alternative, one may resort 
to the diagrammatic theory and to the cooperon-diffuson expansion. A preliminary analysis 
suggests that the expression (83) is of the order E 2 (L/l) 2 (A 2 /jE c ) where L is the circumfer- 
ence and I the mean free path of the metallic ring. Furthermore, 7 the rate at which phase 
coherence is destroyed. Since < Q(fi,y)Q(fi,y') >~ E 2 , the estimate above indicates that in 
metals of short mean free path the stochastic fluctuations may be considerably larger within 
the Fermi level pinning ensemble, that is, when calculated for a fixed number of electrons 
as compared with calculations for fluctuating number of electrons. 



VII. DISCUSSION 

We recall the discussion on the existence and on the size of persistent currents in meso- 
scopic metallic rings threaded by a magnetic flux <fi which took place just a few years ago 
[2]. It has been found that for non-interacting electrons, the persistent currents are expo- 
nentially small [13] when calculated for the grand canonical ensemble (GCE) in contrast 
to calculations for the canonical ensemble (CE). There one finds persistent currents, (when 



23 



expressed in terms of energies) of order A, that is, the mean spacing of the electronic level. 
Nevertheless, even this result is small and it is presumably correct to say that the size of 
the persistent currents depends on just one electron. 

In such a case, it is necessary to have reliable information on the distribution of the 
electronic levels €\. A basic assumption is that atomic disorder or fluctuations in the geo- 
metric definitions exist in all samples which are prepared identically on a mesoscopic scale. 
Correlation functions (with respect to this disorder), of these levels can be calculated in an 
approximation by a diagrammatic method [14] or in the random matrix theory [11,12]. A 
comprehensive approach is provided by the supersymmetry technique [15], but technically, 
it is not so easy to handle. 

Two main subjects are discussed in this paper: (A) Thermodynamics; (B) Dynamics; 
(C) Stochastic Fluctuations. 

Concerning (A), there is no need to emphasize, that for systems of identical particles, 
a theory can be formulated and carried through most elegantly in the grand canonical 
ensemble. Thus, there is a possibility (i) to extract from this information thus obtained the 
canonical ensemble by a Legendre transformation. In a second procedure (ii), one starts from 
the idea of a Coulomb blockade where the charging energy e 2 /2C (C is the capacitance of 
the mesoscopic sample) fixes the number of electrons to be equal to the (effective) number of 
ions provided that e 2 /2C » A. A third possibility (iii) is based on the concept of pinning 
the Fermi level to a single particle level. In this context, one argues that effectively the grand 
canonical ensemble can be understood as a random selection of Fermi energies which prefers 
configurations where the energy separation between last occupied and first empty state is 
large. In contrast to it, the Fermi level pinning ensemble provides an unbiased choice. 

By and large we have found that (at T = 0) the three methods do lead to the same 
result as far as the thermodynamics is concerned. One open question remains. In terms of 
persistent current: Is there a difference in the persistent currents of a flux (p when the system 
has been prepared at a flux 0^0? Some comments can be found at the end of appendix 
A. In this context, the following problem (which at first sight seems to be none at all), 
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should also be considered. We recall Eq. (15) where /x° is defined by —dOP/d/u, ^=^0= N. 
A relation of similar structure is Eq. (34). The question is: how is it possible to find //° 
such that N is an integer? Clearly, within the present formalism, we cannot guarantee N 
to be an integer. On the other hand, such a detail may be very important since we have 
convinced ourselves that it may be only one electron which contributes to the phenomenon 
we are interested in. 

As an example of dynamics (B), we have studied the disorder average cross section 
< a(uo) > for the photoabsorption of a mesoscopic sample as proposed in Ref. [1]. According 
to a qualitative analysis of Shklovskii [1], the cross section < cr(u) > is reduced by a factor 
(ftw/A) (<< 1 for small frequencies) as compared with a naive calculation in the grand 
canonical ensemble. We have calculated < a(uo) > using the concept (iii) of Fermi level 
pinning. There, we have obtained agreement with Shklovskii's ideas within the random 
matrix theory. Furthermore, the present theory allows us to cover all frequencies. 

It seems to be tempting to calculate photoabsorption within the Coulomb blockade con- 
cept (ii). However, in a naive treatment, correlation functions of very high order (— > 00) 
seems to be required; a situation which calls for a more detailed discussion. In this context, 
we wish to mention also an approach to this dynamic problem which has been put forward 
by Kamenev et al. [16], which agrees with the limiting form proposed by Shklovskii [1]. 

Stochastic fluctuations (C) seem to be very important in mesoscopic physics. Though the 
average response may be small, there are in general large sample to sample fluctuations of this 
response. As a simple measure of these fluctuations, the root-mean-square of the response is 
most useful. Usually, the r.m.s is calculated within the grand canonical ensemble. We have 
studied this quantity for Fermi level pinning ensemble (iii). In the case of persistent currents, 
we have outlined an estimate which indicates a considerable enhancement is possible in the 
fluctuations of persistent currents, when calculated for fixed number of electrons. 

At the end, we wish to draw attention to measurements where the magnetic moment 
of (singly connected) mesoscopic samples has been measured by Levy et al. [17]. Again, 
the question of differences between canonical and grand canonical ensemble arises. In the 
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presence of strong disorder (diffusive limit), the paper by Altshuler et al. [18], seems to 
provide a complete answer. On the other hand, if there is only weak atomic disorder, the 
problem is more complicated [19]. A new physical situation arises when samples are free 
from intrinsic defects (Aharonov-Bohm ballistic billiards), and where disorder appears only 
in the form of differences in the geometric definitions of the samples. This situation has been 
discussed by Ullmo et al. [20]. We, together with Yu. N. Ovchinnikov, are also working on 
this problem, analytically as well as numerically. It seems to be a very difficult problem; 
one could say, that it means to catch the one active electron among the million inert ones. 



A.S. wishes to thank A. Kamenev and Y. Gefen for stimulating discussions during his 
stay at the Weizmann Institute. In fact, these discussions have initiated the investigation of 
this paper. In addition, it is a pleasure to acknowledge illuminating comments by A. Mirlin. 
This work has been supported by the German-Israeli Foundation (GIF). 



According to Altshuler and Shklovskii [14], the two level correlation functions can be 
expanded in two types of impurity ladder diagrams that are called diffusons and cooperons. 
For this expansion to be valid, a condition has to be satisfied which is that the energetic 
separation of the two levels to be considered has to be much larger than the mean level 
separation 



In the framework of this diagramatic theory, one obtains a two-diffuson contribution which 
is as follows 
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APPENDIX A: 



E-E' \» A 



(Al) 




(A2) 
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In the relation above, Dcf represents the eigenvalues of the diffusion operator — DV 2 for 
the sample in consideration and 7 means the rate at which phase coherence is destroyed. In 
case of time reversal symmetry, the two-cooperon contributions < T> l {E)T> 1 {E') > 2C is the 
same as expression flA~2f) . 

We call a sample to be effectively of zero dimension if the zero eigenvalue of -DV 2 is 
the only important one. In this case 

<T>\W)>™=-^ {E J l+ . if (A3) 

According to the definition (24), the expression ( |A3| ) has to be considered as the two- 
diffuson approximation of the two- level correlator —y 2 (E,E'). We also observe that this 
approximation satisfies the sum rule (25). Furthermore, it agrees with the GUE result (65) 
of the random matrix theory for large energy separation | E — E' \» A; 7, if we replace 
there the fast oscillatory term sin 2 nx by its average value ~. For small energies, however, 
the agreement is bad. 

As mentioned above, the two-cooperon contribution is the same for a system with perfect 
time reversal symmetry, Y 2 2 c = Y 2D . On the other hand, if time reversal symmetry is 
completely lost, the two-cooperon contribution is zero, Y 2C = 0. In the random matrix 
theory, the two cases correspond to GOE and GUE, respectively. 

Despite its deficiency, the cooperon-diffuson expansion is useful since it can easily be 
generalized to one and higher dimensional samples. In the discussion of persistent currents 
in metallic rings threaded by a magnetic flux </>, one considers frequently an effectively one- 
dimensional closed loop of length L; as a consequence, the eigenvalues Dq 2 of the diffusion 
operator can be found by replacing 

q^-r (n + p) (A4a) 

Lj 

where (p — for diffusons whereas for cooperons 

(p = j- ; 0o = (A4b) 
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We remark that the ansatz (A4) neglects any penetration of the magnetic field into the 
area of the metallic ring. Therefore, the cooperon contribution remains oscillating even for 
very large flux and no transition to the GUE ensemble of random matrix theory will take 
place. 

Presently, we are interested in the thermodynamics only in so far as the flux dependence 
is concerned. Hence, only the two-cooperon contribution Y 2 matters. Considering Eqs. 
(23) and (28), we recognize that the (^-function contribution to Y 2 is irrelevant. We will 
also see later, that neither the 77-regularisation of the low energy dependence nor the phase 
breaking rate 7 is not essential in the present problem. Thus, we may write 

Y 2 2C (u; v?) = At V Re -. -. -—77 7777 (A5) 

2 v ,rj 2tt 2 ^ [u + iD(2ix/L) 2 (n + y?)] 2 V ; 

Note that Y 2 2C is an even function of </?; moreover, it is as periodic in with period 1. 
Therefore, a Fourier expansion is appropriate and we put 

Yi c {u- V )=Y,C m {u)e 2mm * 

m 

/+00 

Explicitly, one finds 



1 1 _ re~^/ 4 .1 m 



c 



Cm{u) = ^^[-^ ~ ^ (A7a) 



where 

E c = D/L 2 (A7b) 

is the Thouless energy and where 

v— u/E c 

a(v)= tt/2 - (7r/4)sgnu (A7c) 
Next, we consider Eq. (28) and calculate the cooperon (if = 20/0 o ) contribution 
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<(5N(<f))) 2 > = J2B m e 



2irim24>/ <fio 



f S 2 1 

B m = s 2 / du uC m (u) = —- 1 (A8) 

Jo ^ I m I 

We observe that the divergence of C m (u) oc u~ 3 / 2 is irrelevant in the integration and also, 

that the integral converges for large values of u if m 7^ 0. This last point may be considered 

as a demonstration of the rule that only low energies (here: | u |< E c ) contribute to phase 

sensitive quantities. 

For sake of completeness, we give the expression for the phase sensitive part of the free 



energy. Inserting expression (|A8|) in Eq. (28), we have 



< AF >= V At— cos 2-nm-r- (A9a) 



m=l 



Since the persistent current 7(0) can be calculated according to 



d 

/ = -c— <AF> (A9b) 
ocp 

we obtain a result that can be expressed as follows 

m=- j{4>) 

JU)= y J?l . 2 ± .^ S in2™^ (A10) 

m=l 

Following the argumentation which lead to Eq. (47), we also define a two-level correlator 
which depends on the two fluxes and 0. In this case, cooperons as well as diffusons 
contribute to the phase sensitive quantities; and we obtain the appropriate correlators if we 
substitute <p = (0 + 0)/0 o for cooperons and <p = (0 — 0) /0 O for diffusons in expression ( |A6|) . 
(Note the symmetry <-» 0.) From these relations, it follows that 



< Afi P (T = 0; 0, 0) > = A < c^V(0)cW(0) > 

f S 2 ^[ + 0, 0-01 / A11 s 

= > — — cos27rm hcos27rm All 



, 7T m 

m=l 



Evidently, 
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/p(0,0)= -c— < AQp(4>,^ > 

= jfy<4> + $!) (A12) 

As expected, the above result agrees with Eq. ( |A10| ) only if = 0. 

Next, we comment the paper of Kamenev and Gefen [10]. They consider an experiment 
where the system is prepared in a grand canonical state at the flux and then transferred 
adiabatically to a state with flux 0. Their result for the persistent current can be expressed 
as follows 

7(0, 0)= -c~A < (6N(<j>) - 5iV(0)) 2 > 

= J{<P) - j(^(0 + 0)) - j(\(<f> - 0)) (A13) 

In conclusion, we note the interesting relation 

Jp(0,0)+ 7(0,0) = J(</>) (A14) 

which we do not understand presently in physical terms. 

As far as the problem of stochastic fluctuations is concerned, it is not difficult to calculate 
the connected part < £7(0)fi(0') > c of the correlator for the grand canonical potential 
in the diffuson-cooperon approximation outlined above (see also the standard results for 
noninteracting electrons quoted in Ref. [3]). Accordingly, one finds for the grand canonical 
ensemble at T = 0, the following result 

< Q 0' > c = V -2-5^ cos cos — (A15) 

^ ™ 5 00 00 

Note that the fluctuations ~ E c are much larger than the mean values ~ A. 

APPENDIX B: 



We outline our calculation of the disorder averaged the normalized cross-section ( |53| ) for 
the GOE. Making use of Eqs. ( |71b| ) and ( |72b| ), we obtain the relation 
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< a > 1 



< a > uc x 



dy Y3(y,y + x,o) 

-~ fdy YM+(l + ^j (l-n(ar)) 
In the RMT, the n point cummulant function is given by the expression 



(Bl) 



:b2) 



where Xij = Xi — Xj and J2 P denotes the sum over the subsets of the symmetric group with 
respect to the group of cycle permutation. We note that in the present problems, these 
permutations produce analytic forms which are identical. 
The matrices cr(x) are given by 

\ 



a(x) 



^ six) D(x) ^ 



s(x) 



dx 



six) 



J(x) s(x) I \ s(t)dt — ^sgn(x) s(x) j 
where s(x) is found in Eq. (|65b|) . For the two point cummulant we get 



(B3) 



Y x (x) = s 2 (x) - D{x)J(x) 

The three point cummulant can be calculated according to 

Ys(y, y + x,0) = 2s(x)s(y + x)s(y) 

+s(x)D(y + x)J{y) + s(x)J(y + x)D{y) 
+D(x)s(y + x)J{y) + D(x)J(y + x)s(y) 
+J(x)D(y + x)s(y) + J(x)s(y + x)D(y). 



(B4) 



(B5) 



We insert these quantities in equation ([Bl]) and after some manipulations, we obtain the 



limiting cases, Eq. (74j) . Numerical calculation have been done accordingly. 
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FIGURES 

FIG. 1. Average cross section < a > of photoabsorption for an ensemble of mesoscopic 
samples, as a function of the photoenergy Hu>, normalized to the mean level distance A (spin 
deceneracy s=2). The cross section is given in the form of the ratio < a > / < a > uc , where 
< o >" c oc to 2 is the uncorrelated quantity. Dotted line: Gaussian orthogonal ensemble; dashed 
line: Gaussian unitary ensemble (GOE and GUE, respectively, of random matrix theory). 
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